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Abstract. It is generally believed that mini black holes decay by emitting elemen¬ 
tary particles with a black body energy spectrum. The original calculation lead to 
the conclusion that about the 90% of the black hole mass is radiated away in the form 
of photons, neutrinos and light leptons, mainly electrons and muons. 

With the advent of String Theory, such a scenario must be updated by including 
new effects coming from the stringy nature of particles and interactions. The main 
modifications with respect to the original picture of black hole evaporation come from 
recent developments of non-perturbative String Theory globally referred to as TeV- 
Scale Gravity. By taking for granted that black holes can be produced in hadronic 
collisions, then their decay must take into account that: (i) we live in a Z?3-Brane 
embedded into an higher dimensional bulk spacetime; (ii) fundamental interactions, 
including gravity, are unified at TeV energy scale. Thus, the formal description of 
the Hawking radiation mechanism has to be extended to the case of more than four 
spacetime dimensions and include the presence of D-branes. This kind of topological 
defect in the bulk spacetime fabric acts as a sort of “ cosmic fly-paper ” trapping 
Electro-weak Standard Model elementary particles in our (3 -I- l)-dimensional uni¬ 
verse. Furthermore, unification of fundamental interactions at an energy scale many 
order of magnitude lower than the Planck energy implies that any kind of fundamen¬ 
tal particle, not only leptons, is expected to be emitted. 

A detailed understanding of the new scenario is instrumental for optimal tuning of 
detectors at future colliders, where, hopefully, this exciting new physics will be tested. 
In this article we review higher dimensional black hole decay, considering not only the 
emission of particles according to Hawking mechanism, but also their near horizon 
QED/QCD interactions. The ultimate motivation is to build up a phenomenologi¬ 
cally reliable scenario, allowing a clear experimental signature of the event. 
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1 Introduction 

Recent developments, both in experimental and theoretical high energy physics, seem to 
“conjure” the astonishing possibility to produce a mini black hole in laboratory. From 
one hand, the next generation of particle accelerators, as LHC, will reach a center of 
mass energy of 0(10) TeV; on the other hand, non-perturbative String Theory, including 
O-branes, can accommodate large extra dimensions characterized by a compactihcation 
scale many order of magnitude bigger than the Planck length. A promising realization 
of these ideas is represented by “TeV-Scale Gravity” where it is possible to introduce 
a unification scale, M*, as low as some TeV (m. i)- If these theoretical models are 
correct new exciting phenomenology (i. i) is going to be observed at future colliders, 
as the physics of unified interactions, including quantum gravity, will become accessible. 
No doubts, the most spectacular expected event is a mini black hole production and 
decay [S]. It has been conjectured ( 0 . 0 ) that two partons colliding with center 
of mass energy \/§ and impact parameter smaller than the effective black hole event 
horizon radius ru{'/I) will gravitationally collapse. One expects that proton-proton 
scattering at LHC will produce, among many possible final states, a rotating, higher 
dimensional, black hole jH] with mass Mbh = ~ 10 TeV ; the estimated rate of 

black hole production in geometrical approximation (i, 01, HH), is as high as one per 
second (1I2], |ISI), less optimistic evaluations give 10^ — 10^ per year. The evaporation 
mechanism considers mini black hole as semi classical objects emitting particles with 
a black body energy spectrum jini 

I ^ P 

{N)^rns = - (- 1 ) 2 * ‘ 

Here, (N )ujms is the average number of particles with energy cn, spin s, third compo¬ 
nent of angular momentum m, emitted by a black hole with angular velocity VLh and 
temperature Tbh- |Ap is the “greybody factor” accounting for the presence of a po¬ 
tential barrier surrounding the black hole. All the relevant quantities (cross sections, 
decay rates, etc.) are compute in the framework of relativistic quantum field theory, 
which provides a low energy effective description of black holes and particle interactions. 
From this point of view. String Theory appears to provide “only” the general scenario, 
including extra dimensions and D-branes, rather than computational techniques. On 
the other hand. String Theory introduces the concept of minimal length as the minimal 
distance which can be physically probed (HEi.ini). Including such feature in the dy¬ 
namics of black hole evaporation leads to a significant change in the late phase of the 
process, as we shall discuss in the conclusion of the paper. 

The paper is organized as follows: in Section |21 we recall the different phases of black 
hole evaporation process, mainly focusing on Schwarzschild black hole decay. In Section 
Elwe consider the QED/QCD interaction among particles directly emitted by black holes 
via Hawking mechanism; thus, in Subsection Id. II we report some numerical results about 

^This is strictly valid only for black holes with masses Mbh much above the fundamental higher 
dimensional Planck scale M* (see El)- 

^Throughout this paper we use the natural units h = c = ks = 
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formation and development of photosphere and chromosphere around a Schwarzschild 
black hole. In Section we analize in more detail parton fragmentation into hadrons 
and black hole emergent spectra: in Subsection 14.11 we consider the case in which a 
chromosphere forms and develops, while in Subsection 14.21 the case in which simple near 
horizon parton hadronization occurs. In Section El we conclude with a brief summary of 
the results and some open problems. 

2 Black Hole Decay 

By taking for granted that a mini black hole can be produced in hadronic collisions, it 
will literally “explode” much before leaving any kind of direct signal in the detectors 

From a theoretical point of view we can distinguish three main phases in the decay 
process: 

1. Spin Down phase: during this early stage the black hole loses most of its angular 
momentum but only a fraction of its mass; numerical simulations Hi indicate that 
75% of angular momentum and about 25% of mass are radiated away. Thus, more 
than one half of the mass is emitted after the black hole has reached a non-rotating 
conhguration. 

2. Schwarzschild phase: during this intermediate stage a spherically symmetric, 
non-rotating, Schwarzschild black hole evaporates via Hawking radiation, losing 
most of its mass. 

3. Planck phase: this is the hnal stage of evaporation, when the residual mass 
approaches the fundamental Planck scale M*, i.e. Mbh — and quantum 
gravity effects cannot be ignored. 

A characteristic feature of higher dimensional models, like “TeV-Scale Gravity”, is that 
a D-dimensional black hole can emit energy and angular momentum both in our (3-1-1)- 
dimensional “brane-universe” {brane emission), and in the D-dimensional “bulk” where 
the brane is embedded {bulk emission) ^2]. A simple estimate suggests that half energy 
is lost in the bulk IZQ: a more detailed calculation shows that black holes radiate mainly 
on the brane, even if the ratio between energy emitted on the brane and in the bulk 
is not much greater than one mm)- As the Schwarzschild phase is considered the 
dominant stage, and only the energy (= particles) emitted on the brane can be directly 
observed, we shall focus only on Schwarzschild black hole brane emission, even if the 
black hole “lives” in D = (4 -|- n) dimensions. In this case, equation reads: 

{N)ujs = ^^/Tbh _ 

where 


Tbh = 


n + 1 
471 rn ’ 


^For an alternative scenario see m 


( 3 ) 
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and rn is the event horizon radins [221 given by 

1 

71 + 1 

(4) 

In order to calcnlate power, P, and flnx, F, emitted on brane we proceed as follows. 
First of all we calcnlate the “greybody factor” | A p in eqnation O- The “greybody 
factor” acconnts for the inflnence of spacetime cnrvatnre on particle motion; a particle 
created near the event horizon mnst cross a “gravitational potential barrier” in order 
to escape to inhnity M- Regarding brane emission, the “greybody factor” can be 
calcnlated as a transmission factor across the potential barrier of the (4 + n)-dimensional 
Schwarzschild metric projected on the brane: 


rH = 




M, 


BH 




71+1 


8r(^) 

71 + 2 


ds^ = — 


dt^ + 

1 - 

/rn W+i' 


\ r / 



V r / 




Once |y4p is known (|2S1, [2^], |2E]), it is possible to define a greybody cross section 1271 
as: 

+ (c..) = ^(2j + l)!|A|y (6) 

where j is the total angnlar momentnm of the particle emitted by black hole. 

The resnlts obtained in 1221 is shown in FignreQ we can see that 


• greybody factor is a fnnction of particle energy (a;) and spin (s); 


• greybody factor depends on the nnmber n of extra dimensions, which lead to 
conclnde that it will be possible to “measnre” the nnmber of extra dimensions at 
fntnre colliders m, m, isni); 

• asymptotically the greybody factor tends to geometric optics cross section 


'^g.o. — 


71 + 3 
2 


2/(n+l) 


71 + 3 
71+1 


H 5 


(7) 


becanse of the hniteness of the gravitational potential barrier, i.e. the potential 
barrier has a maximnm of energy Umax given by: 


^max 


1 In + 1 / n + 3 \ 
rny n + 3 \ 2 ) 


( 8 ) 


When a particle is emitted with an energy u > Umax the black hole behaves as an 
ideal black body with area cr^.o.. 

By integrating eqnation m we can obtain the contribntion of spin s. Standard Model 
particle, to P and F: 

• for a; G ^.Umax]-, we integrate thermal spectrnm (|21) with greybody factor 
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Figure 1: Greybody cross section for brane emission from a (4+n)-dimensional Schwarzschild 
black hole, for with spin 0, ^ and 1 particles (Reprinted figures with permission from [22] . 
Copyright SISSA/ISAS 2003). 


• for cu G [oJraax) oo], we integrate thermal spectrum m with optical geometric cross 
section Cgo.- 

Thus, we find 


P = 



and 


27r2 


dujuj" 


a 


is) 
gb ' 


:a;) 


aio/Tgi 


■1) 


2s 


+ 


27r^ 


duuj^ 


cr, 


g.o. 


olo/Tbi 


-1) 


2s 


( 9 ) 


^ (dN\ 1 , 1 /■“ j 2 

\dt J 2ir^ ij - (-1)^* 2i!^ _ (-1)2* 

( 10 ) 

Accounting for color and spin and considering an equal number of particles and antipar¬ 
ticles emitted by black hole, we find the power and flux emitted by a (4-|-n)-dimensional 
Schwarzschild black hole on brane through different type of Standard Model particles, 
as reported in Table UJ 

We see that ~ 75% of decay products are quarks, anti quarks and gluons, while only 
~ 12% are charged leptons and photons, each particle carrying hundreds GeV of energy. 
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Power 

n = 0 

n = 2 

n = 4 

n = 6 

Quarks 

64,9% 

62,9% 

61,8% 

61,3% 

Charged Leptons 

10,8% 

10,5% 

10,3% 

10,2% 

Neutrinos 

5,4% 

5,2% 

5,2% 

5,1% 

Photons 

1,4% 

1,7% 

1,8% 

1,9% 

Gluons 

11,3% 

13,5% 

14,5% 

14,9% 

Weak bosons 

4,3% 

5,1% 

5,4% 

5,6% 

Higgs 

1,9% 

1,0% 

1,0% 

1,0% 


Flux 

n = 0 

n = 2 

n = 4 

n = 6 

Quarks 

66,5% 

64,1% 

62,1% 

60,8% 

Charged Leptons 

11,1% 

10,7% 

10,4% 

10,1% 

Neutrinos 

5,5% 

5,3% 

5,2% 

5,1% 

Photons 

1,2% 

1,5% 

1,8% 

1,9% 

Gluons 

9,3% 

12,4% 

14,1% 

15,2% 

Weak bosons 

3,5% 

4,6% 

5,3% 

5,7% 

Higgs 

3,1% 

1,3% 

1,2% 

1,2% 


Table 1: Power and Flux emitted by a (4+n)-dimensional Schwarzschild black hole on 
brane through Standard Model particles. 


We conclude that the black hole decay is dominated by partons; since they cannot be 
directly observed, we must to take into account fragmentation into hadrons. Therefore, 
in next section we are going to study the effects of QED/QCD interactions (not only 
hadronization) among particles emitted near event horizon, in order to understand what 
kind of spectra we can expect to detect. 

3 QED and QCD Effects 

In previous section we have outlined the evaporation process of a mini black hole pro¬ 
duced at colliders, with a main focus on the “Schwarzschild Phase”, concluding that 
brane black hole emission is dominated by quarks and gluons. This emission, which 
we shall call direct emission”, consists in the near horizon creation of Standard Model 
particles through Hawking Mechanism and in their flowing to inhnity. However, the pic¬ 
ture consisting of Standard Model elementary particles freely escaping to inhnity with a 
black body energy spectrum must be somehow improved by including QCD interaction 
effects as black hole emission is dominated by quarks and gluons. The inclusion of their 
interactions is instrumental to build up phenomenologically reliable model predicting 
the form of the spectrum to be looked for. 

In the next section we shall consider in more detail parton fragmentation into 
hadrons, here we are going to discuss something which could happen before hadroniza- 
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tion, i.e. the appearance of a quark-gluon plasma around the black hole (|SI1, IS21)- 
One starts by considering quarks and gluons emitted according to a black body law. It 
follows that the parton density near the event horizon grows as If temperature is 
high enough, i.e. above some critical value, one expects from QCD quarks and gluons 
to interact through ^^bremsstrahlung” and pair production” processes. These reactions 
are 2 —> 3 body processes increasing the number of quarks and gluons nearby black 
hole and leading to a kind of quark/gluon plasma surrounding the event horizon. While 
propagating through this plasma quarks and gluons lose energy. When the average en¬ 
ergy is low enough partons fragment into hadrons. 

Similar arguments can be applied to photons, electrons and positrons as well, provided 
we replace QCD interactions with the corresponding QED processes. Thus, we can 
dehne two distinct regions surrounding the black hole as follows: 

• Photosphere, is the spatial region around the black hole where “QED-brems- 
strahlung” and “QED-pair production” lead to formation of an e ^,7 plasma; 

• Chromosphere, is the spatial region around black hole where “QCD-bremsstrahl- 
ung” and “QCD-pair production” lead to formation of a quark/gluon plasma. 

Spectrum distortion induced by the two “atmospheres” dehned above will be discussed 
in the remaining part of this section. 

An analytic description of photosphere and chromosphere dynamics is not available 
at present; the best one can do is to resort to numerical resolution of Boltzmann equation 
for the interacting particle distribution function |221- In order to get a hrst insight of 
the phenomena we are considering, we shall introduce a further simplification: we shall 
use “flat spacetime” diffusion relations even if we are going to study the problem in the 
geometry of a Schwarzschild black hole. 

The number dn{p,r) of particles per unit volume element with momentum between p 
and p + dpis 

dn (p, r) = / (I p I , r) d^pr"^ dr dQ , 

where the distribution function / (|p | ,r) solves the Boltzmann equation 

{P ,Pt C) = ^C[f {p ,pt ,r)] . ( 11 ) 

In fjll|l we introduced the radial velocity Vr, the transverse component of the momentum 
Pi and defined p = |p|. C[f] is the collision term built out of scattering cross sections 
encoding information about interactions taking place inside photo and chromosphere. 

According to QED (QCD), electrons and positrons (quarks and antiquarks) can 
lose energy through the emission of photons (gluons): —»• (gg ^ QQg)- 

Corresponding Feynman diagram is: 
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In the ultra-relativistic limit, differential cross section in the center of mass frame reads 


a 


da 


duj EuirrFj 




X 




V m?jUj 


( 12 ) 


where mj is the fermion mass, a is the electromagnetic (strong) coupling costant, E is 
the initial energy of each fermion, and uj is the energy of emitted photon (gluon). 

To avoid infra-red divergences, it is convenient to use energy-averaged total cross section 


(E5; |Sn]): 


a{E) 



da 

duj 


du! 


~ — log 

nif 



(13) 


The same type of cross section is obtained for pair production, —»• e^e+e {qg 

qqq)] corresponding Feynman diagram is: 



At order one finds: 


cr(o^) ^log "j , (14) 

where ut is the incoming photon (gluon) energy. 

Since a oc 1/m‘j, heavy fermions minimally affect photo/chromosphere development, and 
justify photo/chromosphere description in terms of electrons, positrons and light quarks 
alone, as in pTH] . 

Scattering processes discussed above do not occur in vacuum, rather they take 
place in an hot plasma of almost-radially moving particles, which is what we call 
photo/chromosphere. A simple way to take into account hnite temperature effects con¬ 
sists in replacing vacuum fermion masses (m/ in (j13j) and (HI) with their thermal 
counterparts 

= + , (15) 

where T is the plasma temperature and M is often referred to as the plasma mass. 

The position m gives the propagator pole in momentum space with an accuracy of 
10% [nn]- For finite temperature gauge theories, one finds 

M^{T) = g^C{R)'^, 

where g is gauge coupling constant, C{R) is the quadratic Casimir invariant for gauge 
group representation R. Relevant values of C{R) in numerical computations [221 
C{R) = 1 and C{R) = 4/3 corresponding to U{1) and SU{3) fundamental representa¬ 
tion, respectively. In the hot plasma scenario, fermions move inside photo/chromosphere 








TeV Mini Black Hole Decay at Future Colliders 


8 


as “free” particles with a temperature dependent effective mass mth. 

In next subsection we shall report main results for a numerical resolution of Boltzmann 
equation (HH). 


3.1 Numerical Results 

Let us start this subsection by investigating formation and development of photo/chro¬ 
mosphere around a four-dimensional Schwarzschild black hole [32] , then we shall discuss 
as photo/chromosphere features depend from the number of extra-dimensions. 


3.1.1 Photosphere. 

We already mentioned that scattering processes become important only beyond some 
critical temperature 

Let us introduce Af{r) as the number of collisions a typical particle undergoes between 
the event horizon and some larger radius r > ru- A/'(r) can be written in terms of the 
bremsstrahlung and pair production mean free path, A(r), as 


A/’(r) 



We say that a photosphere surrounds the black hole if every particle is scattered at least 
once between the event horizon and infinity: 


lim A/”(r) > 1. (16) 

r—^oo 

Thus, the critical temperature is the black hole temperature giving one for the 

limit (fT^. From numerical analysis we get 

~ 50 GeV. 


Another important photosphere parameter is the inner radius r* which can be defined 
hj Af (ri) = 1, i.e. the mean radial distance for a particle to be scattered once. Data 
fit (Figure |2|) gives 


n = K = (6.446 ± 0.003) X 10"^ (17) 

The inverse dependence from the temperature is due to the fact that bremstrahlung and 
pair production mean free path decreases with the temperature. Thus, the inner edge 
of the photosphere is closer to the event horizon. By inserting rn = ^t^Tbh m we 
get: 

47r 

rj = —rn — 2 x lOV//. (18) 

K 

Even in the case D > 5, we would not expect a different behavior for inner radius. 
Indeed, bremsstrahlung and pair production are scattering processes involving Standard 


Hn ini, -45.2 GeV. 
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• • • Results 

— Inverse law regression 
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Figure 2: Radii for inner photosphere surface for different temperatures; remember that 1 
GeV“^ ~ 0.197 fm (Reprinted figure with permission from jddj . Copyright 1999 by The 
American Physical Society). 

Model particles bound to the DS-brane. 

From an experimental point of view, it is important to know the average hnal energy Ef 
of particles when they start propagating without signihcant interactions, i.e. the particle 
average energy at outer photosphere surface, because this is the expected energy to be 
eventually released in a detecting device. Figure El shows Ef for different black hole 



Figure 3: Particle average energy at outer photosphere surface versus temperature for Tbh < 
1.5 TeV (Reprinted figure with permission from j33j . Copyright 1999 by The American Physical 
Society). 

temperature Tbh- InTableEl we report some values extrapolated from Figure El Since, 
in the hrst approximation, particles are emitted with a black body spectrum, the average 
energy of a particle close to the event horizon is given by Ei ~ TTbh'-, thus, we notice 
that for Tbh = 60 GeV, Ef is not much different from Ei, because the temperature 
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Tbh (GeV) 

60 

300 

1000 

Ef (GeV) 

97.0 

79.5 

44.2 


Table 2: Average final energy for some values of temperature (Reprinted table with permission 
from m- Copyright 1999 by The American Physical Society). 


is only a little higher than In this regime the rate of bremsstrahlung and pair 

production processes is not high enough to decrease in a signihcant way particle energy. 
On the other hand, for Tbh = 1000 GeV, Ef is much smaller than its initial value. These 
results display in a clear way the role of the photosphere: above a critical temperature 
the number of particles grows and the average energy decreases. This effect is enhanced 
at higher black hole temperature. Figure | 3 ] shows up this behavior in a quite evident 
way. 



Figure 4: Particle number versus energy for a 4D Schwarzschild black hole with Tbh = 1000 
GeV. (a) : near the horizon. (b): at the outer photosphere surface (Reprinted figures with 
permission from m- Copyright 1999 by The American Physical Society). 


3.1.2 Chromosphere. 

By following the same procedure adopted for the photosphere, we hnd 

rpQCD _ ^75 

This result agrees with the analytic estimate in jH^. Thus, we conclude that whenever 
black hole temperature is high enough to produce a photosphere, then a chromosphere 
must be present as well. In such a case, the chromosphere inner radius is close to the 
horizon, i.e. r* ~ ru'- a black hole with temperature Tbh ^ emits interacting 

quarks and gluons in the strong coupling regime with initial black body average energy 
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Ei ~ 3Tbh- As they propagate towards infinity their average energy decreases below 
Aqcd and then fragment into hadrons. In this way, one finds the average final energy 
for quarks and gluons to be 

Ef ~ 200 - 300 MeV. 

Transition from partons into hadrons “marks” the position of chromosphere edge. Chro¬ 
mosphere emerging particle spectrum is fairly different from the direct emission spec¬ 
trum, as there are many more particles with a lower average energy (Figure Ej). 



0.0 3.0 6.0 9.0 12.0 15.0 

Particle energy, GeV 



Figure 5: Particle number versus energy for a 4D Schwarzschild black hole with Tbh = 1-5 
GeV. (a) : Hawking emission spectrum at the horizon and (b) : final spectrum emerging from 
chromosphere (Reprinted figures with permission from | . Copyright 1999 by The American 
Physical Society). 


3.1.3 Concluding Remarks. 

Main results obtained in this section can be summarized as follows: 

• the presence of a photosphere, or a chromosphere, surrounding the event horizon 
implies a proliferation of emitted particles; energy conservation leads to a lower 
average energy per particle. Thus, direct emission spectrum is modified: a black 
hole with horizon temperature Tbh = 1.5 GeV, effectively behaves as black body 
at temperature about 100 MeV (Figure EI). 

When emitted particle interactions are properly accounted for, the “free” Hawking 
spectrum is shifted to a an effective black body spectrum to a temperature lower 
than the black hole temperature. 

• Whenever Tbh ^ Aqcd, black hole emits quarks and gluons and direct emission is 

partons dominated (Table [T]); furthermore, since < Aqcd, strongly coupled 

quarks and gluons form a chromosphere surrounding the event horizon. Thus, one 
concludes that final particle spectra will be dominated by hadrons, coming from 
parton confinement, and their decay products. 
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Results reported in this section have been obtained in the geometry of a four-dimensional 
Schwarzschild black hole. However, examining a more general H-dimensional case, the 
following features have to be taken into account: 

• scattering processes introduced previously are H-independent, because they in¬ 
volve only 3-Brane confined Standard Model particles; so cross sections dn and 
(HI are unaffected; 

• brane emission is an intrinsically four-dimensional process, thus, photo/chromo¬ 
sphere dynamics, emitted power, emitted particle number, average initial energy, 
near horizon particle density don’t change in higher dimensions; 

• brane emission is dominated by quarks and gluons independently of D, further¬ 
more flux percentages of partons and photons, electrons and positrons, are quite 
constant. 

All these considerations lead us to conclude that the four dimensional picture of black 
hole decay does not significantly change when moving to higher dimensions. 

4 Hadronization and Emergent Spectra 

We have seen in previous sections that final spectra, to be measured by an asymptotic 
observer, are going to be dominated by hadrons, coming from parton confinement, and 
their decay products, mainly photons, neutrinos and e^. For this reason, we are going 
to consider first quarks and gluons hadronization and then black hole emergent spectra 
in two cases: i) a chromosphere forms and develops; ii) simple near horizon parton 
hadronization occurs (as in ra). 

In order to recall some usefull formula to study hadronization and hadron decay 
processes, we are going to determine the decay photon emergent spectrum. 

Let us consider the formation of a neutral, light, 7 r° meson. Light mesons, like vr'^, 
represent very likely decay products for heavy hadrons. Furthermore, the preferential 
decay channel, vr^ —>• 7 -|- 7 , produces two photons which can be easily detected. 
Hadronization is an intrinsically non-perturbative effect. As such it is difficult to describe 
in the framework of QCD. Our ignorance about color non-perturbative dynamics can 
be parametrized in terms of the so-called hadronization function. For a neutral pion the 
hadronization function reads j33] : 



(19) 


where gjniQ, L^tt) denotes the number of 7 r° produced in the energy range [Et,, Ej^ + dEj,] 
by an energy Q parton of the j-th kind. Thus, the flux of neutral pions emitted by a 
black hole, per unit time, in the energy range E^^ -|- dF^,r] can be written as 



( 20 ) 
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where the sum runs over all the parton species relevant to neutral pion formation. 

As mentioned above, we can follow two approaches to clarify the physical meaning of 
the term dNj / dQ dt: 

• Direct Hadronization: denotes the number of j-th species partons emitted 

near the black hole horizon, per unit time, in the energy range between Q and 
Q + dQ, according with a black body spectrum 


dN, ^g(Q) 

dQdt TT^ qQ/Tbh _ (^_X)2s ’ 

where s is the parton spin, and is the grey-body cross section (Figure [T]). 
Once emitted near the horizon,quarks and gluons freely propagate toward inhnity. 
However, when the relative distance becomes higher than the threshold value 
^QCD ~ 1 then hadronization starts. 


• Hadronization after chromosphere formation: differently from the preceeding case, 
once emitted near the horizon, quarks and gluons propagate toward infinity form¬ 
ing a chromosphere, as discussed previously. Thus, denotes the flux of the 
j-th species partons near the outer boundary of the chromosphere; this flux can 
be evaluated following the method shown in reference I33|. 


Then, the number of emergent photons, per unit time, in the energy range [E^, E^ + dE^] 
is given by 

dN, r dg,^{E^,EQ dN^ 
dE^dt dE^ dE^^dt ’ 

2 

where, Eq = E^ + is the minimum pion energy which is necessary to produce a 

photon with energy E^, and is the vr® mass. The function number 

of energy E^ photons produced by the decay of a pion with energy Et^ in the laboratory 
frame: 


d^'y-K^^E-iY, EQ 

dE^ 


"7 \/El - m. 

By taking into account equations (Hi, (Ei, (Ei and 


one obtains 


(23) 


dN^ 

dEydt 



dE-n 


15 


SEl^^^El-ml 



dN, _ 

dQdt 


(24) 


This integral can be computed numerically once the partonic flux is given according 
with either of the approaches discussed above. 

Equation ()24|1 provides the total flux of photons which can be experimentally detected; 
these photons are not emitted by the black hole itself but come from the pion decay. 
Indeed we must to them add the near horizon “direct” emission spectrum. However, 
also the direct emission photons can lead to the formation of a photosphere through 
electromagnetic interaction with and produce a modified spectrum (Figure | 3 ]). In 
what follows, we shall list some results obtained for photon emergent spectra. 
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4.1 Hadronization “Post-Chromosphere” 

In this subsection we shall focus on photon emergent spectra accounting hadronization 
after chromosphere formation. 

Figure ini shows the total spectrum of photons emitted by a AD Schwarzschild black hole 
with temperature Tbh = 50 GeV. Both direct emission photons (QED) and 7r° (QCD) 
decay photons are displayed; each case is considered both when photo/chromosphere is 
present (solid curve), and when it is absent (dashed curve). The whole photon spectrum 
is obtained in either case by summing the solid and dashed lines. From Figure El one 



Figure 6: Photon emission spectrum for a 4D Schwarzschild black hole with temperature 
Tbh = 50 GeV. The continuous line denote the spectrum obtained by taking into account 
the presence of photosphere (QED) and chromosphere (QCD). The dashed lines denote either 
the direct emission spectrum, or the (QCD) decay spectrum following near horizon quark 
hadronization (Reprinted figure with permission from ^3]. Copyright 1999 by The American 
Physical Society). 

sees that: 

• peak energy of 7r° decay photons corresponds to an energy equal to m^/2; the 
photon distribution is widened by “Doppler effect” as the pions do not decay at 
rest; 

• according with the previous remark, the photon spectrum in the absence of chro¬ 
mosphere is wider, as the partons do not lose energy in the chromosphere before 
hadronization; 

• the direct emission spectrum is picked around ~ 5Tbh, while photosphere produces 
a larger number of photons with smaller energy; 

• in the QED sector the photon peak is many orders of magnitude smaller than the 
photon peak in the QCD sector; as quarks and gluons dominate direct emission, 
the QCD degrees of freedom, leading to pions, decaying into photons, are many 
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more than direct emission photons. Thus, we conclude that 7 r° decay photons 
dominate the emergent spectrum in Figure 0 both in presence and absence of 
photo/chromosphere. 

Figure [7| shows the decay photon spectrum in the presence of chromosphere (solid line), 
and the direct emission spectrum (dashed line) from a Schwarzschild black hole m D = 
10 dimensions, for M* ~ 1.3 TeV e Mbh/M^, = 5. In this case Tbh — 200 GeV and the 
mean life is tbh — 4x 10“^^s. A qualitative analysis of Figure [7|shows that the two curves 



Figure 7: Chromosphere emergent photon decay spectrum for a Schwarzschild black hole in 
= 10 dimensions, with M* ~ 1.3 TeV and Mbh/M^ = 5. The dashed line shows the 
direct emission spectrum (Reprinted figure with permission from |37j . Copyright 2003 by The 
American Physical Society). 

are peaked around (solid line) and ~ FTbh (dashed line), that is, the involved 

physical processes (direct emission, chromosphere development, hadronization and 7 r° 
decay) are spacetime dimension independent. This is consistent with the hypothesis 
that elementary particle interactions are all localized on the spacetime 3-brane and do 
not feel the presence of bulk extra dimensions. 

We can estimate the number of 7 r° decay photons emitted at the outer boundary of the 
chromosphere to be roughly about 5 ■ 10^ photons E7|. This number can be improved by 
adding the photons emerging from the photosphere and the contributions from different 
hadronic decays. 

4.2 Direct Hadronization 

In this subsection we shall consider some results regarding final energy spectra account¬ 
ing direct hadronization. 

In order to determine emergent spectra of stable species (e^, (w)e/iT and in particular 7 ), 
production and decay of a H-dimensional Schwarzschild black hole has been numerically 
performed by using black hole event generator Charybdis vl.OOl [SH], while both the 
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quark/gluon hadronization and hadronic/leptonic decays have been nnmerically simu¬ 
lated by Pythia v6.227 [32] • 

In Table El we have reported some parameters obtained for 100 generated events. For 


D 

{Mbh) (TeV) 

{Tbh) (GeV) 

(N^) 

(N,) 

(P) 

6 

10.345 (±0.030) 

140.20 (±0.13) 

19.03 (±0.23) 

1483 (±26) 

78.4 (±1.4) 

7 

10.279 (±0.030) 

235.62 (±0.17) 

13.06 (±0.22) 

1244 (±24) 

96.4 (±1.9) 

8 

10.31 (±0.03) 

328.43 (±0.19) 

10.80 (±0.18) 

1138 (±22) 

107.1 (±2.3) 

9 

10.299 (±0.027) 

416.28 (±0.18) 

9.57 (±0.18) 

1028 (±22) 

109.5 (±2.5) 

10 

10.327 (±0.029) 

498.20 (±0.19) 

8.75 (±0.16) 

1014 (±26) 

116.9 (±2.7) 

11 

10.259 (±0.024) 

575.37 (±0.16) 

8.19 (±0.16) 

958 (±22) 

120 (±3) 


Table 3: Direct hadronization parameters obtained by 100 generated events. (Mbh) and 
{Tbh) are the average mass and temperature of a Schwarzschild black hole produced by 
Charybdis generator; {Ni) is the average number of particles directly emitted by black hole 
according to the Charybdis simulation, while (Nf ) is the average number of emergent stable 
particles after Pythia simulated direct hadronization. (P) is the average production factor, 
i.e. the number of emergent stable particles following from each particle directly emitted by 
the black hole. 

different spacetime dimensions D, the simulations show that the average black hole mass 
{Mbh) is approximatively costant, while the average temperature {Tbh) grows; as a 
consequence less and less particles are emitted out of the event horizon but they have 
more and more energy. Thus, the average number of particles directly emitted by black 
hole ( Ni) and the average number of emergent stable particles {Nf) decrease with D, 
while the average number of emergent stable particles following from each particle di¬ 
rectly emitted by the black hole (P) is an increasing function of D. Accounting direct 
hadronization, we can observe that black hole decay has very high multiplicity, i.e. a 
black hole emits a great number of stable particles of O(IO^). 

In Figure IHl we show emergent energy spectra obtained from numerical simulations. 

We notice that the hnal emergent spectra are dominated by photons and neutrinos. Since 
future collider detectors are not tuned to capture neutrino signals, a good signature could 
be missing energy or momentum; thus, in Table E] we show the average missing energy 
{^) and the average missing transverse momentum {pr) for several D. 

Furthermore, in Figure El we show photon energy spectra; since these spectra could be 
directly observed at future colliders, the whole picture provides a possible experimental 
signature of TeV mini black hole evaporation. 

5 Conclusions 

In this article we have reviewed TeV mini black hole decay in models with large extra¬ 
dimensions, by taking into account not only the emission of particles according to the 
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Particle Energy (GeV) 


Figure 8: Stable particle energy spectra for a D = 10-dimensional Schwarzschild black hole 
accounting direct hadronization. 


D 

6 

8 

10 

{^) (TeV) 
{fr> (TeV) 

2.65 (±0.07) 
1.82 (±0.05) 

2.63 (±0.08) 
1.81 (±0.07) 

2.55 (±0.10) 
1.76 (±0.08) 


Table 4: Average missing energy and average missing transverse momentum for a D- 
dimensional Schwarzschild black hole accounting direct hadronization. 

Hawking mechanism (referred to as “direct emission”), but near horizon QCD/QED 
interactions, as well. 

We have focused on higher dimensional Schwarzschild black hole decay and we have ob¬ 
served that “brane emission” is parton dominated (see Table P). Since partons cannot 
be directly observed, we must take into account fragmentation into hadrons; therefore, 
in order to understand what kind of spectra we can expect to detect, we have reported 
emergent photon spectra (Figures El El El and El) , both in the case of more realistic near 
horizon QCD interactions (parton bremsstrahlung/pair production and after that frag¬ 
mentation into hadrons) and in the case of “direct hadronization”. Thus, one hnds that 
hnal emergent spectra, to be measured by an asymptotic observer, are dominated by 
hadrons and their decay products, mainly neutrinos and photons, both in presence and 
absence of photo/chromosphere. 

In the latter case, we have reported some results obtained by using Charybdis/Pythia 
event generator package: black hole decay event is characterized by a large multiplicity, 
as high as 10^ (see Table ED and a large missing energy and missing transverse momen¬ 
tum, as high as TeV (see Table ED- The whole picture provides a possible experimental 
signature at future colliders, as LHC and beyond. 

However, much work has still to be done to obtain a phenomenologically reliable sig- 
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Photon Energy (GeV) 


Figure 9: Photon energy spectra for a H-dimensional Schwarzschild black hole accounting 
direct hadronization. 

nature at collider experiments. For instance, recoil effects of the produced black hole 
im, or the influence of Planck phase on the experimental signature (EH. EH) remain 
still to be accounted for. In the latter case, according to several theoretical frameworks, 
it has been argued that the hnal stage of black hole decay is characterized by a “rem¬ 
nant” formation, i.e. either the black hole temperature abruptly drops to zero EH or 
increases up to a maximum temperature and then continuosly approaches an extremal, 
degenerate conhguration at a hnite black hole mass The effects of this 

remnant formation on black hole evaporation have been investigated in E2|. The main 
result is that black hole emits a larger number of Standard Model particles with a lower 
average energy and transverse momentum than in the case of total evaporation; more 
in detail, the total trasverse momentum {^Pt) is lowered by a quantity of the order of 
the remnant mass. 

In conclusion, the main feature of black hole decay is the missing of energy or transverse 
momentum. When a remnant is left after evaporation the missing transverse momentum 
is of the order of the remnant mass. The big acceptance of the detecting device at LHC 
will enable a complete event reconstruction and to determine the missing energy. A 
precise estimate of the missing energy in the framework of a specihc regular black hole 
model is currently under investigation by the authors and the results will be reported 
elsewhere. 
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